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ABSTRACT 

We present results from the multi-wavelength study of XMMU J1230.3+1339 at z ~ 1. We 
analyze deep multi-band wide-field images from the Large Binocular Telescope, multi-object 
spectroscopy observations from VLT, as well as space-based serendipitous observations, from 
the GALEX and Chandra X-ray observatories. We apply a Bayesian photometric redshift 
code to derive the redshifts using the FUV, NUV and the deep U, B, V, r, i, z data. We make 
further use of spectroscopic data from F0RS2 to calibrate our photometric redshifts, and 
investigate the photometric and spectral properties of the early-type galaxies. We achieve an 
accuracy of Az/(1 + z) — 0.07 (0.04) and the fraction of catastrophic outliers is rj — 13 (0) %, 
when using all (secure) spectroscopic data, respectively. The i - z against z colour-magnitude 
relation of the photo-z members shows a tight red-sequence with a zero point of 0.935 
mag, and slope equal to -0.027. We observe evidence for a truncation at the faint end of the 
red-cluster-sequence and the Butcher-Oemler effect, finding a fraction of blue galaxies fb « 
0.5. Further we conduct a weak lensing analysis of the deep 26'x 26'r-band LBC image. 
The observed shear is fitted with a Single-Isothermal-Sphere and a Navarro-Frenk- White 
model to obtain the velocity dispersion and the model parameters, respectively. Our best 
fit values are, for the velocity dispersion osis = 1308 + 284 km/s, concentration parameter 
c = 4.0^ 4 and scale radius r s = 345*52 kpc. From a 38 ks Chandra X-ray observation we 
obtain an independent estimate of the cluster mass. In addition we create a S/N map for the 
detection of the matter mass distribution of the cluster using the mass-aperture technique. 
We find excellent agreement of the mass concentration identified with weak lensing and the 
X-ray surface brightness. Combining our mass estimates from the kinematic, X-ray and weak 
lensing analyses we obtain a total cluster mass of M™ = (4.56 ± 2.3) x 10 14 M Q . This study 
demonstrates the feasibility of ground based weak lensing measurements of galaxy clusters 
up to z ~ 1 . 



Key words: galaxies: clusters: individual: XMMU J1230. 3+1339 - galaxies: fundamental 
parameters - X-rays: galaxies: clusters - gravitational lensing. 



1 INTRODUCTION 

Clusters of galaxies are receiving an increasing interest in astro- 
nomical studies in the last decade for a variety of reasons. They are 
the largest gravitationally bound objects in the universe and provide 
ideal organized laboratories for the cosmic evolution of baryons, in 
the hot Intra Cluster Medium (ICM) phase and galaxies. They are 
also ideal tracers of the formation of the large-scale structure. 

*E-mail: mike@usm.lmu.de, mlerchst@mpe.mpg.de 



The combination of optical, spectroscopic, lensing and X-ray 
studies allows us to get important insights into physical properties 
of the cluster members, the ICM and the determination of the grav- 
itational mass and its distribution with independent methods. The 
dynamical state of the ICM can also reveal information about the 
recent merger history, thus providing an interesting way to test the 
concept of hierarchical cluster formation, as predicted in the cold 
dark matter (CDM) paradigm. 

In the past, clusters of galaxies have been many time the 
subject of combined weak lensing and dynamical studies (e.g. 
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Radovich et al. 2008), weak and strong lensing plus X-ray stud- 
ies (e.g. Kausch et al. 2007), and weak lensing and X-r ay stu dies 
(e.g. Hoekstra et al. 1 20071) . In additio n, Tay lor et al. 1 120041) for 
ground based data and Heymans et al. (2008) for space based data 
showed for the z = 0.16 supercluster A901/2 (COMBO- 17 survey) 
that a three-dimensional (3D) lensing analysis combined with pho- 
tometric redshifts can reveal the 3D dark matter distribution to high 
precision. 

Wide-field instruments, such as the ESO Wide-Field Imager 
WFI at the ES02.2m with a field of view (FOV) of 34'x33', the 
MegaPrime instrument mounted at the 3.6m CFHT with a FOV 
of «1 square degree, the Suprime-Cam with a FOV of 34'x27'at 
the 8.2m Subaru telescope and more recently the Large Binocular 
Camera (LBC) instruments mounted at the 8.4m LBT (Hill & Sali- 
nari 1998) with a FOV of 26'x26', are particularly well suited for 
the multi-band and weak lensing studies of clusters. One can ob- 
tain multi-band photometry for evolution studies based on spectral 
energy distribution (SED) fitting techniques and one can also apply 
weak lensing studies to the clusters, because the FOV is capable of 
sampling the signal well beyond the cluster size. 

In this paper we present an combined optical, kinematic, weak 
lensing, and X-ray analysis of XMMU J1230.3+1339. This cluster 
of galaxies is located at a = 12h30ml8s, 6 = +13°39'00"(J2000), 
redshift z d = 0.975 and Abell richness^] class 2 (see Fassbender et 
al. 2010, Paper I hereafter, for more details). The cluster was de- 
tected within the framework of the XMM-Newton Distant Cluster 
Project (XDCP) described in Fassbender et al. fc008l) . 

We combine optical 6-band data, multiobject spectroscopy 
data, space based UV data (GALEX) and X-ray (Chandra) data. 

In this work we present a detailed analysis of the optical 
data, estimating photometric redshifts and performing a SED clas- 
sification, exploring the colour-magnitude and colour-colour rela- 
tions. Based on this we conduct a weak lensing analysis and es- 
timate independently and self consistently the cluster mass from 
X-ray, kinematics and weak lensing, complementing the analysis 
presented in Paper I. In section[2]we give an overview of the data 
acquisition and reduction. We briefly outline the reduction of the 
LBC data and describe the acquisition of the spectroscopic redshifts 
in Sec. [2] In Sec.[3]we explain the method used to derive the pho- 
tometric catalogues and redshifts, we discuss the accuracy of the 
derived redshifts for the cluster, and we show the distribution of its 
galaxies. In Sec.|4]we present the weak lensing analysis (mass dis- 
tribution) and constraints for the clusters, whereas the estimates for 
the kinematic analysis (velocity dispersion) are shown in section 
[5] The X-ray analysis (surface-brightness profile, global tempera- 
ture) is presented in section[6] and the mass estimates are shown in 
section [7] We discuss and summarise our findings in section [8] and 

El 

Throughout the paper we assume a ACDM cosmology with 
CIm = 0.27, fl A = 0.73, and Ho = 72 km s Mpc~' . Therefore one 
arcsecond angular distance corresponds to a physical scale of 7.878 
kpc at the cluster redshift. 



2 DATA ACQUISITION AND REDUCTION 
2.1 Optical ground based data 

In Table 1 we summarise the properties (e.g., instrument, filter, ex- 
posure time, seeing) of all optical data used in this analysis. The 
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limiting magnitude quoted in column 8 is defined as the 5-tr detec- 
tion limit in a 1.0" aperture via: 

m, im =ZP- 2.5 * Iog(5 yjN pix (r sk y), (1) 

where ZP is the magnitude zeropoint, N pjx is the number of pixels 
in a circle with radius 1.0" and cr s ^ the sky background noise vari- 
ation. Since this equation assumes Poissonian noise whereas the 
noise of reduced, stacked data is typically correlated, this estimate 
represents an upper limit. 



2.1.1 LBC-observations and Data Reduction 

We observed XMMU J1230.3+1339 in 6 optical broad bands (U, 
B, V, r, i, z), with the 8.4m Large Binocular Telescope (LBT) and 
the LBC instruments in "binocular" mode in an observing run from 
28.0 2 2009 until 02.03.2 009. T he LBC cameras (see e.g. Giallongo 
et al. l2008l , Speziali et al. l2008l) have a 4 CCD array (2048 x 4096 
pixel in each CCD; 0.224"pixel scale; 26'x26' total field-of-view). 
Three chips are aligned parallel to each other while the fourth is 
rotated by 90 degrees and located above them. The U,B,V images 
were obtained with the blue sensitive instrument (3500 - 9000 A) 
and in parallel to the r,i,z images with the red sensitive instru- 
ment (5000 - 10000 A). The filter transmission curves are shown 
in Figure [TO] For the single frame reduction and the creation of 
image stacks for the scientific exploitation, we follow the method 
described in detail in Rovilos et al. ( 120091) . Here we quickly sum- 
marise the most important reduction steps: 

1. Single frame standard reduction (e.g. bias correction, flat- 
fielding, defringing). 

2. Quality control of all 4 chips of each exposure to identify bad 
ones (e.g. chips with a large fraction of saturated pixels, etc.). 

3. Identification of satellite tracks by visual inspection and cre- 
ation of masks. 

4. Cr eation of weight images with Weightwatcher (Bertin & 
Marmo 120071) for each chip, accounting for bad pixels and the 
satellite track masks. 

5. Ex tractio n of catalogues with SExtractor (Bertin & 
Arnouts [19961) for astrometric calibration. 

6. Astromet ric c alibration with the USNO-A2 astrometric cata- 
logue (Monet Il998l) using the IRAF ccmap package. This yields 
an external accuracy with an rms of 0.25 arcsec a nd an internal 
accuracy with an rms of 0.066 arcsec (Rovilos et al. 2009). 

7. Flux calibration with standard star observations in each of our 
6 bands using Landolt standards or SDSS-DR7. 

8. Background subtraction and measure of the variation of cr^y 
of each exposure to identify bad images (e.g. those with a value 3 
cr above the median value). 

9. PSF control of all exposures to identify bad patterns (e.g. 
strong variation of the direction and amplitude of the measured star 
ellipticity on each chip of each exposure, etc.) 

10. Co-addition with the Swarp software (Bertin 1200 8l) . 

11. Masking of the image stacks to leave out low-SNR regions 
and regions affected by stellar diffraction spikes and stellar light 
halos as well as asteroid tracks. 

The number of frames which remain in this procedure is dif- 
ferent from band to band, see column 5 in Table 1. Starting from 
50,33,20,20,33,50 images in the U,B,V,r,i,z filters, we at the end 
only stack 24,17,13,12,27,44 frames for the photometric and weak 
lensing analyses. In Figure [4] we show as an example one "good" 
and one "bad" single frame exposure. 
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Figure 1. XMMU J1230. 3+1339 at z = 0.975: 6'x6'LBT-i band image (5.94 ks), overlaid with colour coded contours. In yellow: the surface brightness 
contours of the smoothed light of cluster members, selected to lie at z c luster - 0.05 < z p hot °S ^-duster + 0.05; in blue: weak lensing M a( , significance contours 
(1 - 5); in red: X-ray contours (0.3 - 2.4 keV); in magenta: VLA-FIRST (0.3 - 0.9 mJy) contours. The red big box marks the location of the BCG, the red 
small boxes mark the high confidence spectroscopic confirmed members (13), the blue small boxes (numbers are redshift and error) mark the photo-z selected 
cluster members, and the black circles mark a radius of 0.5 and 1 Mpc from the cluster centre, respectively (note: r2oo~ 1 Mpc). The large object to the west is 
the foreground galaxy NGC 4477. 



Table 1. Cluster optical data overview. 



Field/ Area 
[sq. deg.] 


Telescope/Instrument 


FOV 
[arcmin*arcmin] 


Filter 


N" 


expos, time 
[s] 


zero-point 
[AB mag] 


miim * 
[AB mag] 


seeing 
["] 


C"sky 

[ADU] 


XMMUJ1230.3+1339 


LBT/LBC B 


26*26 


U 


24 


4320 


32.22 


26.97 


0.85 


3.2 


(0.19) 


LBT/ LBC B 


26*26 


B 


17 


3060 


32.98 


26.90 


0.87 


6.9 




LBT/LBC B 


26*26 


V 


13 


2340 


33.27 


26.32 


0.85 


15.3 




LBT/LBC R 


26*26 


f 


12 


2160 


33.44 


26.18 


0.95 


20.3 




LBT/LBC R 


26*26 


i 


27 


4860 


33.27 


25.80 


0.87 


24.5 




LBT/LBC R 


26*26 


z 


44 


7960 


32.79 


25.37 


0.78 


23.5 



Note: a number of single exposures (t cxp = 180 s) used for the final stack; * the limiting magnitude is defined in Eq.[T] c lensing band. 
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Figure 3. The 6'x6', three-colour image of the cluster XMMU J1230. 3+1339, with r-channel (5.94 ks i-band), g-channel (3.24 ks r-band) and b-channel 
(5.94ks B-band) colours. The cluster members appear as reddish objects. The large, yellowish object on the right side is the foreground galaxy NGC 4477. 
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Figure 4. LBT i-band single frames: the left panel shows a "good" frame 
(FWHM 0.64") whereas the right panel shows a "bad" one (FWHM 1.20"). 

We derive our number counts by running SExtractor on the 
final image stacks with the corresponding weight map. We con- 
volve the data with a Gaussian kernel and require for a detection 
to have 4 contiguous pixels exceeding a S/N of 1.25 on the con- 
volved frame. Stars are then identified in the stellar locus, plotting 
the magnitude versus the half-light radius in Figure[20]and identify 
high signal to noise but unsaturated stars. In Figure[5]we show the 
achieved number counts for the final image stacks. Table [T] sum- 
marises the properties of each image stack. For comparison we plot 
in Figure [5] fo r the U , B, V band, the number counts derived from 
Rovilos et al. for the deep Lockman Hole observation with 

the LBT. We do not show the uncertainties for these literature re- 
sults, which have large errors at the bright end (U, B, V mag < 
22.5). In addition we compare our B, V, r and i band number counts 
to those from McCracken et al. |2003) for the VVDS-Deep Field. 
For the z band t he comparison is with the number counts from Er- 
ben et al. J2009I) for the CFHTLS Dl Field. In terms of amplitude 
and slope, our agreement with the literature measurements is very 
good. The comparison of our data to deeper literature data and the 
linear number-count fits added with a solid line shows that the 50% 
completeness for an extented galaxy roughly falls together with the 
prediction for a point source within a 2" diameter aperture (long 
dashed line) according to Equation 1 . 

2.2 Spectroscopy 

In Table 2 we give an overview of all spectroscopic data used in 
this analysis. 

2.2.1 VLT Spectroscopy and Data Reduction 

For the kinematic study of the cluster and the calibration of the 
photometric redshifts we use the data taken in April 2008 and June 
2008 with the FOcal Reducer and low dispersion Spectrograph 
(F0RS) mounted at the UT1 of the Very Large Telescope (VLT), 
which were taken as part of the VLT-program 081.A-0312 (PI: 
H. Quintana). The MOS/MXU sp ectrosc opy was performed with 
FORS2 (see e.g. Appenzeller et al. 1 19981) in medium resolution (R 
= 660) mode (grism GRIS_300I+21). The data were reduced using 
the IRArQ software package. For the redshift determination (red- 
shift, error, confidence level) the spectra were cross-correlated with 
template spectra using the IRAF package RVSAO. More details can 
be found in Paper I. 

The spectroscopic observations provide 38 spectra, where 17 
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of them matched in position with objects from our photometric cat- 
alogue. For the comparison to photometric redshifts we consider 
galaxies with trustworthy > 95% confidence spectroscopic redshifts 
only. Due to low S/N (at high redshift) and the limited wavelength 
ranges of the spectra only 8 objects could be considered. 

We have been awarded additional FORS2 spectroscopy for ~ 
100 photo-z selected cluster members; the analysis of this data will 
be described in an upcoming paper. 



2.2.2 SDSS Spectroscopy 

The field of XMMU J 1230.3+ 1339 has complete overlap with the 
Sloan Digital Sky Surve^ (SDSS). Via the flexible web-interface 
SkySer veiQ of the Catalog Archive Server (CAS), we get access 
to the redshift catalogue. For our purpose we only consider objects 
clearly classified as a galaxy with a redshift > 95% confidence. By 
matching the SDSS catalogue with our catalogue we find 13 objects 
with spectroscopic redshifts between and 0.6 in the SDSS-DR7 
(Adelman-McCarthy et al. 2007). For the comparison to photomet- 
ric redshifts we consider 10 of the 13 objects for which the confi- 
dence level is larger than 3. 



2.3 GALEX UV data 

We include near-UV (NUV, 1770-2730 A) and far-UV (FUV, 1350- 
1780 A) data from the Galaxy Evolution Explore^] (GALEX). The 
instrument provides simultaneous co-aligned FUV and NUV im- 
ages with spatial resolution 4.3 and 5.3 arcseconds respectively, 
with FOV of a; 1.2°. The zero-points have been taken from GALEX 
Observer's GuidH ZP FUV = 18.82 and ZP NUV = 20.08 

mag AB . 

We use data taken during the nearby galaxy survey 
(NGS) and the guest investigator program (Gil) and pub- 
licly released in the GR4. In total we make use of three 
pointings, GI1_109010.NGC4459, NGA_Virgo_MOS09, and 
NG A.Virgo JVIOS 1 1 , with an total exposure time of « 8.4 ks in 
NUV and ss 3.0 FUV. After retrieving the data from the web- 
interface MAST0, we further process them with public available 
tool<0 (SCAMP, SWarp) and obtain resampled and stacked images. 

The GALEX data are rather shallow. Nevertheless they are 
useful to improve the quality of the photometric redshifts, in par- 
ticular to disentangle SED degeneracies in the blue part of the spec- 
tra as detailed in Gabasch et al. (2008). Figure [6] shows that some 
galaxies can clearly be detected in the GALEX filters (FUV,NUV), 
e.g. foreground spiral galaxy above the brightest cluster galaxy 
(BCG) shows an NUV-excess. 



3 Funding for the SDSS and SDSS-II has been provided by the Alfred P. 
Sloan Fundation, the Participating Institutions, the National Science Foun- 
dation, the U.S. Department of Energy, the National Aeronautics and Space 
Administration, the Japanese Monbukagakusho, the Max Planck Society, 
and the Higher Education Funding Council for England. The SDSS Web 
Site is http : //www . sdss . org/ 

4 http : //cas . sdss . org/astro/en/tools/search/SQS . asp 

5 http://galex.caltech.edu 

6 http: //galexgi . gsfc.nasa.gov/docs/galex/Documents/ERO_data_description_; 

7 http://galex.stsci.edu/GR4/ 

8 http://astromatic.iap.fr/ 
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Figure 5. Number counts per 0.5 mag and sq. degree of galaxies (crosses) and stars (open triangles) for the image stacks in the 6 (U, B, V, r, i, z) filters. 
The black line shows a linear fit to the galaxies number counts. For comparison we plot the number counts (open circles) derived from Rovilos et al. (2009) 
and (filled circles) McCracken et al. (2003), see text for details. The short/long dashed lines indicate the estimated limiting magnitudes respectively, defined 
as a point source 5<x detection limit in a 1.0/2.0" aperture. We used Eq.[T]for the noise levels of the image stacks listed in the Table 1. By comparing our 
numbercounts to deeper literature data and the extrapolation of the linear fit we can estimate the 50% completeness for the detection. We have shaded the 
region where the completeness drops below 50% in each filter. 




2.4 X-ray data 

2.4.1 Chandra observations 

We use the archival observation by Chandra (PI: Sarazin, Obs-Id 
9527), as listed in Table 3. The cluster centre lies on the S3 chip 
of the Advanced CCD Imaging Spectrometer (ACIS) array. The 
observation was conducted in VFAINT mode, and reprocessed ac- 
cordingly using the Chandra softwar^l CIAO 4. 1 (Friday, Decem- 



ber 5, 2008) and related calibration files, starting from the level 1 
event file. We used the CIAO 4.1 tools ACIS_PROCESS.EVENTS, 
ACIS_RUN_HOTPIX to flag and remove bad X-ray events, which 
are mostly due to cosmic rays. A charge transfer inefficiency (CTI) 
correction is applied during the creation of the new level 2 event 
filS This kind of procedure reduces the ACIS particle background 
significantly compared to the standard grade selectiorF*! whereas 



9 http://cxc.harvard.edu/ciao 



id 
1 1 



http : //cxc .harvard . edu/ciao/threads/createL2/ 

http : //asc . harvard . edu/cal/Links/Acis/acis/Cal_prods/vfbkgrnd 
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Table 2. Spectroscopic data overview. 



ID 


RA 


DEC 


Mag_Auto 


7 ° 


h 

Flag£, 


ref/_ 


comments 




faegj 




lAB 










34 


187.5598939 


13.7739387 


17.1514 


0.1300 


4 


1 




35 


187.6228679 


13.7854572 


16.8121 


0.1573 


4 


1 




44 


187.5528796 


13.4550020 


16.5786 


0.0830 


4 


1 




46 


187.5733437 


13.4931093 


16.3311 


0.0832 


4 


1 




47 


187.5878466 


13.4371208 


17.0668 


0.0838 


4 


1 




50 


187.6041016 


13.4623009 


16.2996 


0.0829 


4 


1 




57 


187.3367446 


13.6621547 


20.2306 


0.3918 


4 






58 


187.4747995 


13.6681400 


17.4096 


0.1447 


4 






59 


187.4367670 


13.5817725 


18.3299 


0.4377 


4 






61 


187.6479182 


13.6215654 


18.3519 


0.2439 


4 






67 


187.5400466 


13.6773148 


24.3561 


0.9642 


4 


2 


[Oil] emitter 


69 


187.5498491 


13.6759957 


22.0181 


0.9691 


4 


2 




80 


187.5740730 


13.6573756 


21.9994 


0.9758 


4 


2 




82 


187.5681276 


13.6472714 


20.8572 


0.9786 


3 


2 


BCG 


83 


187.5767988 


13.6519471 


22.2796 


0.9793 


3 


2 




84 


187.5791935 


13.6507665 


22.2532 


0.8616 


3 


2 




85 


187.5838292 


13.6503877 


22.8933 


0.9690 


4 


2 




86 


187.5851446 


13.6441070 


21.9350 


0.9767 


4 


2 





a Measurement error: <t z ~ 0.0002. 

" Flags and confidence classes - 0: no redshift, 1: 50 % confidence, 2: 75 
% confidence, 3: 95 % confidence, 4: 100 % confidence 
c 1: Adelman-McCarthy et al. 12007); 2: Paper I 



Table 3. Cluster X-ray data overview. 



Object 


Telescope 


Instrument 


Filter/Mode total exposure clean exposure 

[s] 


Net counts" 


Obs-Id 


comments 


XMMU J123S. 3+1339 


Chandra 


ACIS-S 


VFAINT 37681 36145 


639 


9527 


chip S3 



Notes: " number of net detected counts in the 0.5 - 2 keV band. 



source X-ray photons are practically unaffected (only about 2% of 
them are rejected, independent of the energy band, provided there 
is no pileup 12 !). Upon obtaining the data, aspect solutions were ex- 
amined for irregularities and none were detected. Charged particle 
background contamination was reduced by excising time intervals 
during which the background count rate exceeded the average back- 
ground rate by more than 20%, using the CIAO 4.1 LC.CLEAN 
routine, which is also used on the ACIS background file;0. Bad 
pixels were removed from the resulting event files, and they were 
also filtered to include only standard grades 0, 2, 3, 4, and 6. After 
cleaning, the available exposure time was 36145 seconds. 

Images, instrument maps and exposure maps were created in 
the 0.5 - 7.0 keV band using the CIAO 4.1 tools DMCOPY, MKIN- 
STMAP and MKEXPMAP, binning the data by a factor of 2. Data 
with energies above 7.0 keV and below 0.5 keV were excluded due 
to background contamination and uncertainties in the ACIS cali- 
bration, respectively. Flux images were created by dividing the im- 
age by the exposure map, and initial point source detection was 
performed by running the CIAO 4.1 tools WTRANSFORM and 
WRECON on the 0.4-5.0 keV image. These tools identify point 

12 http : //cxc . harvard . edu/ciao/ahelp/acis_pileup . html 

13 http : //cxc . harvard . edu/ciao/threads/acisbackground/ 



sources using "Mexican Hat" wavelet functions of varying scales. 
An adaptively smoothed flux image was created with CSMOOTH. 
This smoothed image was used to determine the X-ray centroid of 
the cluster, which lies at an RA, Dec of 12:30: 17.02, +13:39:00.94 
(J2000). The uncertainty associated with this position can be ap- 
proximated by the smoothing scale applied to this central emis- 
sion (2.5"). Determined in this manner, the position of the X-ray 
centroid is 15± 2.5 " (110.85+ 18.5 h^ 1 kpc) from the BCG of 
XMMU J1230.3+1339. Figure [2] shows an LBT image of XMMU 
J 1230.3-t- 1339 with adaptively smoothed X-ray flux contours over- 
layed. 



2.5 Radio data 

For the identification of X-ray point sources we make use of Faint 
Images of t he Ra dio Sky at 20-cm from the FIRST survey (see 
Becker et al. ll995h . Using the NRAO Very Large Array (VLA) and 
an automated mapping pipeline, images with 1.8"pixels, a typical 
rms of 0.15 mjy, and a resolution of 5"are produced, and we re- 
trieve this data through the web-interfacj^l 



http : //third . ucllnl . org/cgi-bin/f irstcutout 
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Figure 7. Photometric error as a function of apparent magnitude in the 6 
(U, B, V, r, i, z) filters. 

Figures □ and |2] show LBT images of XMMU J1230.3+1339 
with smoothed radio contours overlaid, see Paper I for more details. 



3 OPTICAL ANALYSIS 

In this section we give an overview of the catalogue creation, the 
photometric redshift estimation, and the galaxy classification. Fur- 
thermore we describe in detail the construction of luminosity maps 
and the colour-magnitude relation. 



3.1 Photometric Catalogues 

For the creation of the multi-colour catalogues we first cut all im- 
ages in the filter U, B, V, r, i, z to the same size. After this we mea- 
sure the seeing in each band to find the band with the worst seeing. 
We then degrade the other bands with a Gaussian filter to the see- 
ing according to the worst band. For the object detection we use 
SExtractor in dual-image mode, using the unconvolved i-band 
image as the detection image, and measure the fluxes for the colour 
estimations and the photometric redshifts on the convolved images. 
We use a Kron-like aperture (MAG_AUTO) for the total magnitudes of 
the objects, and a fixed aperture (MAG.APER) of x 2" diameter (cor- 
responding to a 15 kpc at the cluster redshift) for the colours of the 
objects. Figure[7]shows the median photometric error as a function 
of apparent magnitude for the different filters. The magnitude er- 
ror is from SExtractor which does not consider correlated noise. 
Comparing the formal magnitude error with the 50% detection lim- 
its suggests a scaling of a factor of two. To account for Galactic ex- 
tinction we apply the correction factor E(B - V) = 0.022 retrieved 
from the NASA Extragalactic Databas^fl from the dust extinction 
maps of Schlegel et al. ( 1998), to our photometric catalogue. 

To account for extended halos of very bright stars, diffraction 
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Figure 8. Semi-automatic image masking for the 26'x26'LBT field. The 
polygon squares mark objects from the density analysis and the stars cover 
sources identified in the GSC-1, GSC-2 and USNO-A2 Standard Star Cata- 
logues, see Erben et al. (2009) for further details. 

spikes of stars, areas around large and extended galaxies and vari- 
ous kinds of image reflections we generate a mask, shown in Figure 
[8] We do not consider objects within these masks in the photometric 
and shape catalogue later on. 

3.2 Photometric Redshifts 

We use the Bayesian PHOTO-z code (Bender et al. 2001) to estimate 
the photometric redshifts. The method is similar to that described 
by Bemtez et al.(2000). The redshift probability for each SED type 
is obtained from the likelihood according to the x 1 of the photom- 
etry of the object relative to the SED model, and from a prior de- 
scribing the probability for the existence of an object with the given 
SED type and the derived absolute magnitude as a function of red- 
shift. The code has been successfully applied on the MUNI CS sur- 
ve y (Dro ry et al. l2004) . the FORS Deep Field ( Gabasch et al. l2004al 
& I2006I) . the COSM OS (Ga basch et al. l2008l) . the Chandra Deep 
Field (Gabasch et al. l2004bl) . and the CFHTLS (Brimioulle et al. 
120081) . 

Small uncertainties in the photometric zero-point and an im- 
perfect knowledge of the filter transmission curve (see Figure ITOt 
have to be corrected. We use the colours of stars from the Pickles 
(1998) library for these adjustments. We compare the colours for 
this library of stars with the observed colours and apply zeropoint 
shifts to achieve a matching. We fix the i-band filter and shift the 
remaining ones by -0.12,0.02,0.01,-0.15 and 0.10 for the U,B,V,r 
and z band filters. With additional spectroscopic information we 
can improve the calibration for the photometric redshift code by 
slightly adjusting the zero-point offsets for each band again. The 
match of our observed colours after zero-point shift to the Pickles 
(1998) stellar library colours is shown in Figure|9] 

Using all position matchable spectra (31), leads to an accu- 
racy of <x AJ /(i+j) = 0.074 and the fraction of catastrophic out- 
liers is ~ 13%, as plotted in Figure [12] The dotted lines are for 
Zpimi = Zspec ± 0.15 (1 + Zs P ec)- The fraction of catastrophic out- 



The Massive Galaxy Cluster XMMU J 1230.3+1339 atz~ 1 9 




Figure 9. Colour-colour diagrams of stars (black dots), plotted against the Pickles (1998) stellar library (red dots). The stars used in this plot are selected 
with 17.0 < mag < 23.0 in all bands. The photometric errors (as shown in Figure 7) are of the order of 0.02 mag. All colours perform well. To obtain these 
colour-colour diagrams we have applied zeropoint corrections to all the bands, in detail -0.12, 0.02, 0.01, -0.15, 0.00 and 0.10 magnitudes for U,B,V, r, i, and 
z bands. 
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Figure 10. Transmission curves of the 6 (Bessel-like U,B,V from LBCB; 
Sloan-like r,i,z from LBCR) filters. The dotted lines show the CCD effi- 
ciency of the LBCB (left) and LBCR (right). LBCB has the more blue sen- 
sitive curve. The dashed line shows the atmospheric absorption at the ob- 
servatory site. The red lines correspond to the effective filter transmission 
curves, corrected for CCD efficiency and atmospheric absorption. 

liers is defined as rj = \z spec - z pto; |/(l + z spec ) > 0.15. cr Az/(l+z) is 
the redshift accuracy measured with the normalised median abso- 
lute deviation, 1.48 x median(\Az\/(l + z)). Using only the secure 
spectra (17) which fulfill the quality criteria Flag > 3 (> 95 % con- 
fidence, see Paper I) the accuracy becomes cr AZ / ( i + -) = 0.042 with 
no catastrophic outliers. 

In Figure [T3l the galaxy redshift histogram of all objects in 
the field is shown. The galaxy redshif t distri bution can be parame- 
terised, following Van Waerbeke et al. 2001: 

where (zo,a,/3) are free parameters. The best fitting values for the 
cluster field are zo = 0.77, a = 1.0, and [} = 1.6. For the fit we have 
replaced the peak at z = 1 by the average value of the lower and 
higher bin around the cluster. The fit is shown in Figure [T3l 

Since the number of spectroscopic data for this field is still 




phS * 




Figure 11. Redshift probability distribution (upper panel) and photomet- 
ric redshift fit (lower panel) of the spectroscopic object 82, the BCG with 
Zspec = 0.9786. The different SED's are colour coded, the open circles 
represent the SED integrated within the various filter transmission curves, 
whereas the black dots represent the measured fluxes. The redshift for the 
best fitting SED is Zpkot = 0-94 ± 0.04, which is in agreement with the 
spectroscopic one. 

limited we use another method to estimate the photometric red- 
shift accuracy: we measure the two-point angular-redshift correla- 
tion of galaxies in photometric redshift slices with 14.0< i^B <26.0 
selected from the total sample as shown in Figure [T3] and show 
the results in Figure [14] In this way the null hypothesi s of n on- 
overlapping redshift slices can be tested (see Erben et al. 2009 and 
Benjamin et al. 12010] for a detailed description of this technique). 
For the calculation of the two-point correlation function we use 
the publically available code ATHEN/0- The blue curves in Fig- 
ure [14] show the autocorrelation of each redshift-bin. The red data 
points show the cross correlation of different redshift slices, which 
should be compatible with zero in absence of photometric red- 
shift errors. The redshift binning starts from and increases with 
a stepsize of 0.1 up to redshift 1.0, and stepsize of 0.5 from red- 
shift 1.0 to redshift 4. Nearly all cross-correlation functions of non- 
neighbouring photo-z slices show indeed an amplitude consistent 
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Figure 12. Photometric redshifts, plotted against the spectroscopic ones for the cluster field around XMMU J1230. 3+1339. To distinguish the galaxy spectral 
types (see SED template plotted in the right panel) we use a colour-code: red dots represent early type galaxies (e.g., ellipticals, E/SO), yellow, magenta and 
green dots represent spiral like galaxy types (e.g., Sa/Sb,Sbc/Sc), and blue dots represent strong star-forming galaxies (e.g., Sm/Im & SBm). The median 
residual is (z p hot ~ Z s pec)/(1 + Zspec) = 0.042, o" az /(i +z ) = 0.074 and r\ = 0% for all secure matchable spectra of the cluster XMMU J1230. 3+1339. Note that 
the filled circle mark secure spectroscopic redshifts. 



-i — i — | — ■ — i — i — i — r 

Cluster 



t 1 r— 




Figure 13. Redshift number distribution of all galaxies. The number distri- 
bution is shown as black histogram for the 26' X 26' field, with a median 
redshift Zmedian — 0.77. The solid red curve represents a fit to the galaxy 
distribution, using the approach of Van Warbeke et al. (2001). The arrow 
indicates the galaxy cluster. The shaded distribution show the objects used 
for weak lensing analysis lateron. These are all objects for which photo-z 
and shape could be measured within a FOV 26' X 26' and outside the masked 
area of Figure 8. 



with or very close to zero. This is a strong argument for the robust- 
ness of the photometric redshifts. No excessive overlap between 
low- and high-z slices is observed. The weak signal in the high red- 
shift slice is due to the poor statistics in this and the neighbouring 



bins. Note that the last two bins contain less than 200 galaxies each. 
For the weak lensing analysis (Sect.|4j the vast majority of galaxies 
studied is limited to z = [1.0 - 2.5]. The correlation of all redshift 
intervals with galaxies in z e [1.0 - 2.5] is very small. 

3.3 Light distribution 

To quantify the light distribution we compute the absolute K-band 
luminosity for all the galaxies. Based on the SED fitting, performed 
during the photo-z estimation, we estimate for a given filter curve 
(K-band) the flux and the absolute magnitude. Next we map the 
flux of all galaxies of the cluster in a redshift interval Zduster ± 0.05 
on a pixel grid, and smooth it with a circular 2' Gaussian filter. 

Figure[T]shows an LBT image of XMMU J1230.3+1339 with 
surface brightness contours overlayed. The cluster seems to be very 
massive, rich and has a lot of substructure, see also Paper I. 



3.4 Colour-Magnitude Relation 

The colour-magnitude relation (CMR, Visvanathan & Sandage 
1 1977b is a fundamental scaling relation used to investigate the evo- 
lution of galaxy populations. The CMR of local galaxy clusters 
shows a red- sequence (RS, Gladders & Yee 2000) formed of mas- 
sive red elliptical galaxies undergoing passive evolution. For the 
following discussion we consider as cluster member those galaxies 
in a redshift interval Zciusier ± 0.05, within a box of 6'x6'centered 
on the X-ray centroid of the cluster. Our sample is neither spectro- 
scopic nor does it use the statistical background subtraction as it is 
done when purely photometric sample are studied. As shown in De 
Lucia et al. d2007l) for the ESO Distant Cluster Survey (EDisCS), 
results are not significantly affected whether one does use spec- 
troscopic or photometric redshifts or does not apply photometric 
redshifts for member selection. Note that the accuracy of our pho- 
tometric redshifts is better by a factor of 1.5 than those used for 
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Figure 14. Two-point angular-redshift correlation of galaxies in the field of the cluster XMMU J1230. 3+1339 with 14.0 < iab < 26.0 selected from the total 
sample as shown in Figure 13. The inlay shows as an example the auto-correlation for the redshift bin 0.9 < Zphot < 1-0; the scale of the angular distance 
is logarithmic from 0.1 to 10 arcmin, whereas the correlation function scales linearly from -0.25 to 0.85 . The blue curves show the auto-correlation of each 
redshift-bin. The redshift bins start from and increase by a stepsize of 0. 1 until redshift 1 .0 and stepsize of 0.5 from redshift 1 .0 to redshift 4. The red curves 
show the cross-correlation of galaxy subsample pairs at different redshifts. 



EDisCS (see details in Pello et al. l2009l and White et al. l2005l) . This 
is due to the fact that we can use deep 8 band data (FUV, NUV, U, 
B, V, r, i, z) in contrast to 5 band data (V, R, I, J, Ks) for EDisCS. 

The CMR for XMMU J1230.3+1339 is shown in Figure [B] 
To distinguish the cluster galaxy spectral types we use a colour- 
code: red dots represent early type galaxies (e.g., ellipticals, E/S0), 
yellow and green dots represent spiral like galaxy types (e.g., 



Sa/Sb,Sbc/Sc), and blue dots represent strongly star-forming galax- 
ies (e.g., Sm/lm and SBm). 

The CMR shows a tight red sequence of early type galax- 
ies (ETG). The expected colour for composite stellar population 
(CSP) models with old age and short e-folding time scales or of 
old single stellar populations (SS P) is i ^g - Zab ~ 1-0 (e.g., see 
Santos et al. 120091 : Demarco et al. l2007t) . We apply a linear fit to 
the 9 spectroscopic confirmed RS members at zab < 24 and obtain 
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Figure 15. Colour-magnitude relation, i'ab -Zab against zab. FOV 6'x6', the dashed lines show the 5<r limiting magnitudes and colours. The cluster members, 
are selected to lie at Zciuster ~ 

0.05 < Zphot < Zchmer + 0.05. We use the same colour coding as shown in Figure 12. The black solid line refers to the linear 
fit to the passive cluster members, reproduced with red dots. The dotted lines correspond to the 3ct region. The shaded region on the right indicates > 
50 % incompleteness for a point source within a 2" aperture, which is roughly equal the 50% completeness for an extended galaxy. Black boxes mark 
spectroscopically confirmed cluster galaxies. The dashed-dotted limited region marks the area where we select background galaxies for the WL analysis. 



i- z = (-0.078 ± 0.03 1) x (z- 24) + 0.849 ± 0.065, a result with poor 
statistics due to the limited number of spectroscopic data. Therefore 
we also perform a fit to the 49 photo-z RS members at zab < 24 and 
obtain i-z= (-0.0273 ± 0.012) X (z - 24) + 0.935 ± 0.015. We 
take the best-fitting RS for the photo-z RS members and plot it in 
Figure [15] The slope is consistent in the errors with that obtained 
by Santos et al. J2009h for XMMU J 1229+015 1, the "twin brother" 
cluster at the same redshift observed with ACS onboard the HST. 

We can observe spiral and star-forming galaxies (yellow, 
green, blue data points) until zab - 25.5. The ETG (red data points) 
can only by detected until zab = 24. Note that the 50% complete- 
ness in z-band is 24.5 mag and in i-band is 25 mag, which is the 
detection band. We observe 7 galaxies with 23. 4< z < 23.8. If the 
slope of the luminosity function (LF) is a a -1 we expect about 10 
galaxies with 23. 8< z < 24.2. We observe only one in this inter- 
val. To explore whether this lack of red galaxies is only apparent, 



caused by incompleteness of red the galaxy detection, we have plot- 
ted < 50% completeness region for the i and z detection as a shaded 
area in Figure [15] We estimate the detection efficiency of galaxies 
with 23. 8< z < 24.2 and i — z ~ 1 to be « 50%. Hence we expect 
to detect about 5 galaxies against the only one found with 23. 8 < z 
< 24.2. Furthermore we plot the numbercounts of all cluster galax- 
ies and the ETG ones only as a function of z-band and absolute 
V-band magnitudes in Figure [18] We can see that ETG's are less 
and less starting at z = 22.5 mag, for which the data are 100% 
complete. A similar trend can be seen in the absolute V-band mag- 
nitude counts, where ETG's with My < -23.5 become rare. Similar 
results have been reported for the EDisCS clusters in De Lucia et 
al. (2007 ) and the cluster RDCSJ0910+54 at z = 1.1 (Tanaka et 
al.|2008j). Demarco et al. (2007) report an indication for truncation 
of the i775 - zgso versus zsso CMR (see Figure 15 in Demarco et 
al. 2007) at a magnitude of z 8 50 = 24.5 for RDCSJ1252.9+2927 at 
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Figure 16. Colour-colour diagram with J/ab _ #AB against tab - zab- The same FOV, photometric redshift selection and colour coding as in Figure 15 are 
used. On the right panel we show the colour-errors for the early type galaxies. Overplotted we show single stellar population (SSP) models for passively 
evolving galaxies with Salpeter-like IMF and for three different metallicities (1,2 and 3 times solar). The tracks show the time evolution from 100 Myr (filled 
circles) to 6 Gyr in steps of 100 Myr for two redshifts (0.0 and 0.975). Black boxes in both figures mark spectroscopically confirmed cluster galaxies. 
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Figure 18. Number counts for the galaxies within a 6'x6' region centered on XMMU J1230. 3+1339 and a photometric redshift of Zcluster ~ 0.05 < Zphot < 
Zduster + 0.05. The open and shaded histograms show all cluster galaxies and early type cluster galaxies, respectively. Error bars represent the Poisson errors. 
The red short (long-dashed) line represent the 100 (50)% completeness in the i-filter, the detection band, and for comparison in black the z-filter. The dotted 
line represents m(z) ~ 23, at which galaxies are classified into 'luminous' and 'faint'. There is a rarity of red galaxies with z > 22.2 or i > 23.2 or My ^ -24.2 
where the detection efficiency is still 100% in both bands. 



z = 1.24. A simila r trend at zgso-24.0 can be seen in Figure 8 of 
Santos et al. j2009t) for XMMU J1229+0151 at z = 0.975. 



3.5 Blue galaxy fraction - Butcher & Oemler effect 

The fraction of blue galaxy types increases from locally a few per- 
cent, 25 % (z = 0.4) to 40 % (z ~ 0.8), see e.g. Schade Jl997l) . This 
is called the Butcher & Oemler effect. 

Butcher & Oemler jl978h defined as blue galaxies those which 
are by at least A = 0.2 mag bluer in the B - V rest-frame (i - z 
observer-frame) than early-type galaxies in the cluster. They then 
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Figure 17. Absolute V-band magnitudes of galaxies in the observed 
26'x26'field as a function of redshift. The black box shows the region for 
the selection of cluster members later on. 
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Figure 19. The Butcher & Oemler effect in the colour-magnitude relation, 
>ab ~ Zab against My. The dashed lines show the selection box used by 
Butcher & Oemler (1984) - see text for details. We use the same region, red- 
shift selection and colour coding as in Figure 15. The contours are colour 
coded: black lines represent all photometric redshift selected cluster mem- 
bers, red and blue contours represent the early and late type cluster popula- 
tion, respectively. 



The dashed lines in Figure [19] mark the selection boundaries if one 
would follow exactly the work of Butcher & Oemler. The colour 
coding of the galaxy types is shown in Figure [19] We can infer 
from our photometric analysis a fraction of blue galaxies ft = 0.514 
+ 0.13 for the cluster and f b = 0.545 ± 0.04 for the field at the same 
redshift, with Poisson errors. Note that for the estimate of the frac- 
tion of blue galaxies in the field we take the full 26' x 26'FOV into 
account and exclude the region around the cluster (6'x 6')- 
Our result is in good agreement with the prediction from literature. 
This can be seen if one extrapolates the blue fraction estima ted fo r 
the EDisDS sample from z ~ 0.8 to z = 1 in De Lucia et al. d2007h . 
Within the framework of a state-of-the-art semianalytic galaxy for- 
mation model Menci et al. ( 2008) derived an equal fraction of blue 
galaxies for the cluster and the field at z~ 1 . 

3.6 Colour-Colour Diagram 

We now turn our attention to the colour-colour diagram of the 
galaxies. In Figure[l6]we show r—z against U—B colours and over- 
laid single stellar population (SSP) models0f° r passively evolving 
galaxies with Salpeter-like IMF and for three different metallicities 
(1, 2 and 3 times solar). We expect a Uab~Bab colour of « 1.5 from 
old SSP models with solar metallicity, see tracks in Figure[T6] The 
colour-colour diagram shows Uab - Bab values from -0.5 to 2.5 
in the distribution of the elliptical galaxies (red objects in Figure 
1 1 6b - Note that the Uab - Bab colours are affected by large U-band 
measurements errors (see Figure [7]l. Nevertheless, we estimate the 
weighted mean of all ETG in the cluster, which results in a mean 
colour of Uab - Bab = 0.81 ± 0.5 . This could be a mild indication 
for a recent star formation (SF) episode. 

The observed U-B color of a galaxy at z~l corresponds to a 
rest-frame far-UV color. It may be no surprise the fact that early- 
type galaxies at this redshift exhibit a rest-frame far-UV color about 
0.8 mag bluer than expected from an SSP model. In fact, at least 
30% of the early type galaxies in the local Universe seem to have 
formed 1 %— 3% of their stellar mass within the previous Gyr (Kavi- 
raj et al. 120071) . Hence, it is plausible that a larger fraction of the 
stellar mass of an early type galaxy at z~ 1 is formed within a look- 
back time of 1 Gyr. This would produce significantly bluer far-UV 
colors in the rest-frame if attenuation by internal dust is negligible. 
Alternatively, such bluer colors may spot the presence of a signifi- 
cant UV upturn phenomenon (see Yi 2008, for a review). We will 
investigate the nature of the bluer rest-frame UV colors in early- 
type galaxies in clusters at z< 1 in a future paper. 



4 WEAK LENSING ANALYSIS 

In this section we give an overview on the applied shape measure- 
ment technique and on the lensing mass estimators; we use standard 
lensing notation. For a broad overview to the topic, see e.g. Mel- 
lier dl999h . Bartelmann & Schneider d200lh and Hoekstra & Jain 
d2008t> . 



counted red galaxies down to an absolute magnitude M v = -19.3 
mag within an arbitrary reference radius and applied a statistical 
subtra ction of the background. In contrast to Butcher & Oemler 
dl984l) we use a redshift/SED selection as defined in the previ- 
ous section. Note that it has been shown by De Lucia et al. (2007) 
that a photometric redshift selection does not affect the result. As 
seen in Figure [19] the cluster shows a bimodal galaxy population. 



4.1 Gravitational shear estimate with the KSB method 

The ellipticity of a background galaxy provides an estimate of the 
trace-free components of the tidal field cPifr/dOiddj, where tf/(9) is 
the two-dimensional potential generated by the surface density k{0) 

17 Models computed using PEGASE2 (Fioc and Rocca-Volmerange 1997). 
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(Kaiser & Squires 1 19931) of the mass distribution of the cluster. Any 
weak lensing measurement relies on the accurate measurement of 
the average distortion of the shape of background galaxies. There- 
fore it is crucial to correct for systematic effects, in particular for 
asymmetric contributions to the point spread function (PSF) from 
the telescope and from the atmosphere. 

The most common technique is the so called KSB method pro- 
posed by Kaiser et al. ( 1995) and Luppino & Kaiser dl997l) . In the 
KSB approach, for each source the following quantities are com- 
puted from the second brightness moments 



eU,2), 



(3) 



where W r is a circular Gaussian weight function with filter scale 
r g , and 1(6) is the surface brightness distribution of an image: 
(1) the observed ellipticity 
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(2) the smear polarizability P sm and (3) the shear polarizability P sh . 
The KSB method relies on the assumption that the PSF can be de- 
scribed as the sum of an isotropic component which circularises the 
images (seeing) and an anisotropic part introducing a sys temati c 
contribution (for a detailed discussion see e.g. Kaiser et al. 1 19951 . 
Luppino & Kaiser 1 19971 and Hoekstra et al. 119981 - KSB+). 

The source ellipticit y e s of a galaxy after PSF smearing 
(Bartelmann & Schneider 1 200 lh is related to the observed one, e, 
and to the reduced shear g by: 



(5) 



where the anisotropy kernel p describes the effect of the PSF 
anisotropy (an asterisk indicates that parameters are derived from 
the measurement of stars): 



Pi 



(6) 



The p quantity changes with the position across the image, there- 
fore it is necessary to fit it by, e.g. a third order polynomial, so 
that it can be extrapolated to the pos ition o f any galaxy. The term 
P s , introduced by Luppino & Kaiser dl997h as the pre-seeing shear 
polarizability, describes the effect of the seeing and is defined as: 



pS _ psm I V m \ 
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(7) 



with the shear and smear polarisability tensors P sh and P m , calcu- 
lated from higher-order brightness moments as detailed in Hoekstra 
et al. dl998l) . In Figure ?? The computed p^'/p"" values for differ- 
ent bins of r h and r AB are shown in Figure [23] 

Following Hoekstra et al. dl998h . the quantity p sl "/p s »" 
should be computed with the same weight function used for the 
galaxy to be corrected. The weights on ellipticities are computed 
according to Hoekstra et al. (2000), 



(p g r 



(P^a-l + <A e 2 > 



(8) 



where cr^ ~ 0.3 is the scatter due to intrinsic ellipticity of the 
galaxies, and (Ae 2 ) 1 ' 2 is the uncertainty on the measured elliptic- 
ity. In Figure [2J] we show the galaxy weight as a function of the 
absolute first and second components of the ellipticity. 
We define the anisotropy corrected ellipticity as 



We choose r g to be equal to the SExtractor parameter FLUX_RADIUS. 
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Figure 20. The r-band magnitude versus half-light radius. The red zone 
shows unsaturated stars selected for the PSF anisotropy correction (19 < 
Tab < 22.5 mag, 2.1 < rj, < 2.6 pixel), where one pixel equals 0.224 
arcsecond. 
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Figure 21. Galaxy weight (see Eq.[8) as a function of the absolute first and 
second components of the ellipticity for the weak lensing measurement. 



er = e,-P^p), 

and the fully corrected ellipticity e' so as 



(9) 



(10) 



which is an unbiased estimator for the reduced gravitatonal shear 

r 



= g = gi + ig2 
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assuming a random orientation of the intrinsic ellipticity. For the 
weak distortions measured, k <k 1 , and hence, 
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Figure 22. Distribution of the r;, and tab magnitude values of all galaxies used for the weak lensing analysis. The shaded distribution shows the 785 objects 
used to perform the weak lensing study on XMMU J1230.3+1339. 
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Figure 23. PSF correction: P,;, / P„„ values of all galaxies used for the weak lensing analysis computed in different bins of ry, and r^g magnitude and for each 
chip. 
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Our KSB+ implementation uses a modified version of Nick 
Kaiser's IMCAT tools, adapted from the 'TS' pipeline (see Hey- 
mans et al. 2006 and Schrabback et al. 2007), which is based on 
the Erben et al. ( 2001) code. In particular, our implementation in- 
terpolates between pixel positions for the calculation of Q t j, P sm , 
and P sh and measures all stellar quantities needed for the correction 
of the galaxy ellipticities as function of the filter scale r g following 
Hoekstra et al. {l99it). This pipeline has been tested with image 



simulations of the STEP projectfj . In the analysis of the first set of 
image simulations (STEP1) a significant bias was identified (Hey- 
mans et al. 2006), which is usually corrected by the introduction of 
a shear calibration factor 



<g> = c ca i ■ (e iso ) 



(13) 



with c cal = 1/0.91. In the analysis of the second set of STEP image 
simulations (STEP2), which takes realistic ground-based PSFs and 
galaxy morphology into account, it could be seen that with this 



19 http : //www . physics . ubc . ca/$\sim$heymans/step . html 
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Figure 24. PSF anisotropy correction: ellipticities of the stars in the field of XMMU J1230. 3+1339 as observed (top-left panel) and fitted (top-right panel), 
plus their residuals (bottom-left panel) after correction for PSF anisotropies as a function of their pixel position of the r-band stack. The direction of the sticks 
coincides with the major axis of the PSF, the lengthscale is indicating the degree of anisotropy. An ellipticity of 10% is shown in the upper left corner of two 
upper and the bottom left panel. The definition of ellipticity is given in Eq. 4. The stars initially had systematic ellipticities up to ~ 3 - 8% in one direction. 
The PSF correction (bottom-right panel) reduced these effects to typically < 1.5%. 



calibration factor the meth od is on average accurate to the ~ 3% 
leveE] (Massey et al. 120071) . 

Note that for the weak lensing analysis we only consider the r- 
band stack. We do neither use the i-band stack, because we observe 
larger residuals of the anisotropy correction, nor we do make use of 
the z-band stack, which shows remaining fringe patterns after data 
reduction. We select stars in the range 19 mag < r AB < 22.5 mag, 
0.48"< r/, < 0.58", giving 322 stars (see Figure l20l to derive the 
quantities for the PSF correction. Figure [24] shows the distortion 



Note: Dependencies of shear on size and magnitude will lead to different 
uncertainties. 



fields of stellar ellipticities before and after PSF anisotropy correc- 
tion as a function of their pixel position of the r-band stack. The 
direction of the sticks coincides with the major axis of the PSF, 
the lengthscale is indicating the degree of anisotropy. The stars ini- 
tially had systematic ellipticities up to ~ 3 - 8% in one direction. 
The PSF correction reduced these effects to typically < 1.5%. The 
achieved PSF correction for LBC is at the same level in terms of 
the mean and standard deviation of stellar ellipticities as the one 
for Suprime-Cam as detailed in Table 4 and Figure 2 of Okabe & 
Umetsu (2008). 

Galaxies used for the shear measurement were selected 
through the following criteria: w > 3, SNR >5, z pllol > 1.03, r;, > 
0.58", tab < 25.5 mag, and ellipticities, e iso , smaller than one. 
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This yields a total number of 785 background galaxies in the cat- 
alogue, and an effective number density of background galaxies 
(Zphoi > 1-03) of n e ff ~ 7 galaxies per square arcmin. In Fig- 
ure Q3] we show the distributions in the in the i AB - zab against 
Zab colour-magnitude diagram of the members of the galaxy clus- 
ter XMMU J1230.3+1339 and the background galaxies. The area 
is limited with the dashed-dotted line. Note that without accurate 
photometric redshifts a separation of background from foreground 
and cluster member galaxies cannot be achieved. 

The measurement of the shape of an individual galaxy pro- 
vides only a noisy estimate of the weak gravitational lensing signal. 
We therefore average the measurements of a large number of galax- 
ies as a function of cluster centre distance. Figure [26] shows the re- 
sulting average tangential distortion e, of the smear and anisotropy 
corrected galaxy ellipticity as a function of the distance from the 
X-ray centroid of the galaxy cluster. The errorsbars on e, and e c 
indicate Poissonian errors. We detect a significant tangential align- 
ment of galaxies out to a distance of ~ 6'or ~ 2.8 Mpc from the 
cluster centre. The black points represent the signal e c , when the 
phase of the distortion is increased by n/4. If the signal observed 
is due to gravitational lensing, e c should vanish, as observed. To 
obtain a mass estimate of the cluster, we fit a singular isothermal 
sphere (SIS) density profile to the observed shear. We use photo- 
metric redshifts for the computation of angular distances, where 
D s , D d , and D ds are the distances from the observer to the sources, 
from the observer to the deflecting lens, and from the lens to the 
source. We show our results in Figure [26] where the solid line is 
the best fit singular isothermal sphere model which yields an equiv- 
alent line of sight velocity dispersion of a = 1308 ± 284 km s~ ! and 
X 2 IDOF = 1.66/5. 

In addition to the SIS model considered so far from now on 
we use in addition an NFW model (Navarro et al. l 1997T) to describe 
the density profile of the cluster. This is because a NFW model au- 
tomatically provides a virial mass M 2 oo- An NFW model at redshift 
z is defined as 



P(r) 



6 C Pc 



(r/r s )(l+r/r s f' 
where 

c 200 c 3 

Oc 



3 ln(l +c)-c/(l +c) 



(14) 



(15) 



is the characteristic density, c = r 2 oo/ r s is the concentration param- 
eter, r s is the scale radius, and p c (z) = 3H(z) 2 /SnG is the cosmo- 
logical critical density. The mass within a radius R is given by 



M(r<R) = 4n6 cPc (z)rl[ln(l + ^)- T ^-] 
and the corresponding mass Mioo is given by 
M 200 = 200p r (z)^r 2 3 00 . 



R/r s 



(16) 



(17) 



We perform a likelihood analysis (see Schneider et al. 2000) 
for details of the method) to fit an NFW model to the measure- 
ments. This yields a concentration parameter of c = 4+.J 4 , a scale 
radius of r, = kpc and ^ 2 /DOF= 1784/783 within a 50% sig- 

nificance interval, as shown in Figuresl27l 

We test our lensing signal for contamination by the imperfect 
PSF anisotropy correction. Following Bacon et al. J2003h we com- 
pute the correlation between the corrected galaxy and uncorrected 
stellar ellipticities. To assess the amplitude we normalise the quan- 
tity by the star-star uncorrected ellipticity correlation. 




r [ercTnln] 

Figure 25. PSF systematic residual: the cross-correlation function f™ s (red 
points) between galaxies and stars is shown as a function of angular scale. 
The amplitude of the cross-correlation is always at least one order of mag- 
nitude smaller than the shear signal £e ~ e} (black filled). 



fair) : 



(e'(x)y(x + r)) 1 



(18) 



e] up 



(e*(x)e*(x + r)) ' 

where the "asterisk" indicates a stellar quantity. In Figure [ 
show the cross-correlation compared to the shear signal £ £ ~ 
to 10'. Note that the amplitude is at least one order of magnitude 
smaller than the signal. We conclude that the PSF correction is well 
controlled. 



4.2 Aperture mass M ap 

We make use of the aperture mass statistics, introduced by Schnei- 
der d 19961) and worked out further by Schneider d 19971) to quantify 
the mass density concentration and the signal of the weak lensing 
detection. The aperture mass statistic has the advantage that it can 
be derived from the shear field in a finite region and M ap is not 
influenced by the mass-sheet degeneracy. It is defined as, 



d 2 9 k(0) U(0 - 8Q, 



(19) 



where U (6) is a radially symmetric weight function, using a com- 
pensated filter function of radius 0L, 



f 

Jo 



dO 6 U(8) = 0. 



One can express M ap in terms of the tangential shear e, 

M ap (e c ) = J d 1 6e,(9) Gdfl-flgD- 
The filter functions Q and U are related through 
2 



Q{e) = 



e 1 



f d& 9' 
Jo 



U(ff) - U(9). 



(20) 



(21) 



(22) 



It was shown by Schneider dl996h and Bartelmann & Schneider 
j200ll) that the variance cr Map of M ap and M ap itself can be esti- 
mated for observational data in a straightforward way. In fact one 
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Figure 26. Radial profiles of the tangential shear e ( (in red, measure of the 
signal) and the cross shear e c (in black, measure of the systematics) around 
the cluster XMMU J1230.3+1339. The error bars indicate the statistical 
error in the measurement, due to the intrinsic ellipticity of the individual 
sources. The black solid line shows the fit of a SIS model to the data, which 
yields a velocity dispersion of <x = 1308 ± 284 km/s (x 2 /DOF = 1.66/5), 
while the dashed line shows the NFW model. 




Figure 27. Likelihood confidence contours for the NFW fit to the two shear 
components on a grid centered on the cluster. The solid line contours are 
for one, two, and three <x confidence levels (68.3%, 95.4%, and 99.73% re- 
spectively) with best fit parameters c = 4.0* ^ and noo = 1376^2?g kpc 



(X 2 /DOF 
level. 



1784/783). The dotted line indicates the 50% significance 



only has to change the integral into a sum over galaxy ellipticities 
and the error analysis is simple expressed by 



M„ 



" el gj e,j Wj 0m - ggP 

* 2 ^2, 4 W ? 2 2 (ie,' 



2(2, w,y 



(23) 



(24) 



where e, , are the tangential components of the ellipticities of lensed 
galaxies relative to each of the grid points investigated, and w, are 
the weights defined in Eq[8] 

Therefore the signal to noise ratio (S/N) can by calculated as 



S/N-- 



M„ 



(25) 



We use the polynomial filter function introduced by Schneider 
et al. dl996l) : 

6 (em) 2 , , N 
Q(e) = --^{i-mf), (26) 

where 6 is the projected angular distance on the sky from the aper- 
ture centre, and (f Q is the filter radius. We apply this weight function 
to the data. 

We detect a weak lensing signal at a 3.5 cr significance. The lo- 
cation of the highest significance coincides with the X-ray measure- 
ment of the Chandra analysis (see section[6]l. In Figure[T]we show 
an LBT image of XMMU J1230.3+1339 with M ap S/N-contours 
overlayed. We observe an elongation in the projected mass distri- 
bution from weak lensing in the same direction as the distribution of 
galaxies associated with the main cluster component and the group 
bullet (see Paper I), and the X-ray emission. 



5 KINEMATIC ANALYSIS 

For the cluster kinematic analysis the redshifts and positions of 
galaxies are used to define a characteristic velocity dispersion and 
a characteristic size (scale radius) of the cluster. According to the 
virial theorem, the RMS velocity dispersion has the property to be 
completel y indep endent of any shape and viewing angle (see Carl- 
berg et al. 1 1996 



5.1 RMS velocity dispersion 

Following Carlberg et al. 1 1997b the line-of-sight velocity disper- 
sion is defined as 

, Z, Wj vl 



Z,w, ' 

where the w, are magnitude dependent geometric weights, 
_ cfa - id) 



(27) 



(28) 



the peculiar velocity in the frame of the cluster, and z L -i is the mean 

redshift of the cluster. 

An alternative approach (Carlberg et al. [l997) is the un- 
weighted line-of-sight velocity dispersion, defined as 



7 v 2 



N- 1 



(29) 



For a virial mass estimator the projected mean harmonic point- 
wise separation can be used 
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(30) 



where ij indicate a sum over all pairs. The three-dimensional virial 
radius is defined as 



r v tr 



71 R H 



The virial mass of the cluster is defined as 

3 cr; r rir 
M vir = ' . 



(31) 



(32) 



We use the model of a singular, isotropic, isothermal sphere 
with a of velocity dispersion <t S!S = cr,,, and therefore a density 
profile of 



P(r) 



(33) 



IjxGr 2 ' 

To compare masses usually an extrapolation to a constant mean 
interior density of 200 p c „r is performed, using 



Ma 



M w > 



p(r vir ) 



I 200 p c „, 

where the mean density p(r,.„ ) inside r rjr is defined as 

3 M V1> 



P(r vir ) 



4 7r r 3 ,. 



(34) 



(35) 



r | p " x-iv I ' 
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Figure 28. Radial surface brightness profile of XMMU J1230.3+1339 for 
the 0.29-7.0 keV band sampled in 4 arcsec annular bins (pixel size = 0.492 
arcsec/pixel). The red solid line traces the best fitting single fi model with 
r c = 39.99 pixel (155 h~\ kpc) and p = 0.67. 



5.2 Cluster Membership 

Next we define the cluster membership of galaxies, commonly this 
is determined in two ways: 



5.2.1 3o~ clipping technique 

The s tatistical clipping technique (e.g., Zabludoff, Huchra & Geller 
1 19901) is one method to identify cluster members. It is based on the 
measured recession velocities of the galaxies. Yahil & Vidal fl977» 
showed that this is a reasonable assumption, if the velocity distri- 
bution of the galaxies follows a Gaussian distribution and clusters 
are relaxed and isothermal. 



5.2.2 Mass model method 

Another technique to assess the boundaries of a cluster under usage 
of the spatial and redshift data is the mass model method (Carlberg, 
Yee and Ellingson 1 19971) . We make use of the model described 
in Pimbblet et al. ( 120061) for the Las Campanas Observatory and 
Anglo-Australian Telescope Rich Cluster Survey (LARCS), which 
is based on Carlberg et al. (1996). 

For our kinematic mass estimate we use the line-of-sight ve- 
locity dispersion cr™ = 657 ± 277 km s -1 , compiled in Paper I. 

Based on the assumption that a cluster is a singular isothermal 
sphere (Carlberg et al. 1 19971) . 



'200 



10 H{z) 



(36) 



where cr,. ss V3 cr^ is the velocity dispersion of the clus- 
ter and H(z) 2 = Y§[Cl M (\ + zf + O t (l + zf + Sl A ], We estimate 
r 20a = 944 + 397 kpc and M 20 o = (2.85 + 3.60) x 10 14 M Q ; note 
that due to the poor statistics we only give upper limits. We would 
need more spectroscopic data of cluster members to improve the 
estimate together with its uncertainty. 



6 X-RAY ANALYSIS 

6.1 Surface brightness profile 

A radial surface brightness profile was computed from a combined, 
point-source removed image and exposure map over the range 0.29- 
7.0 keV in 4" annular bins, using the CIAO 4.1 tools DMEX- 
TRACT and DMCALC. 

The surface brightness profile was th en fitted with a single /? 
model (Cavaliere & Fusco-Femiano [l976l) . of the form: 



/(r) = / 1 + - 



.2 \ -3y3+ 4 



(37) 



where Iq is the normalization, r c is the core radius, and can be 
determined from the profile at r » r c . 

The best fitting ft model of this cluster is described by r c = 
39.99 ± 12.58 pixel (155 ± 46.5 K£ kpc), /3 = 0.67 ± 0.20 (see 
Figure l28t. The reduced x 2 of the fit was 1.20 for 27 degrees of 
freedom. 



6.2 Global temperature 



Using the X-ray centroid position indicated in Section 12.4.11 a 
point source removed spectrum was extracted from the observa- 
tion in a circular region with a 40" (3 15 h^, 1 kpc) radius. The spec- 
trum was then analysed in XSPEC 11.3.2 (Arnaud 1 19961) . using 
weighted response matrices (RMF) and effective area files (ARF) 
generated with the CIAO 4.1 tools SPECEXTRACT, MKACIS- 
RMF, MKWARF and the Chandra calibration database version 4.1. 
The background was extracted from the aimpoint chip in regions 
free from the cluster emission. The resulting spectrum, after back- 
ground subtraction, extracted from the observation contains 639 
source counts. The spectrum is fitted with a singl e temperature 
spect ral model (MEKAL, see Mewe et al ll985lll985l : Kaastra et al. 
1 19921 : Liedahl et al. 1 19951) including Galactic absorption, as shown 
in Figure lAll The absorbing column was fixed at its measured value 
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Table 4. Comparison of different X-ray statistics. 



Cash-star* 


Cash-star 3 


^ 2 -stat 


^ 2 -stat 


^ 2 -stat 


^ 2 -stat 


(minimally binned, 


(minimally binned, 


(20 cts/bin, 


(20 cts/bin, 


(30 cts/bin, 


(30 cts/bin, 


N/y,Z/Z Q ,z c / fixed) 


N/y,Z c ; fixed) 


N/7,Z/Z ,z c ; fixed) 


N/y,z c ; fixed) 


N/y,Z/Z ,z c ; fixed) 


N/y,z c ; fixed) 





2.64 


2.64 


2.64 


2.64 


2.64 


2.64 


kT 


5.29 ± 1.02 


5.28 ± 1.02 


4.73 ± 1.12 


4.85 ± 1.12 


4.75 ± 1.16 


4.76 ± 1.17 


z/z G 


0.30 


0.29 ± 0.24 


0.3 


0.74 ± 0.43 


0.30 


0.54 ± 0.43 


^cluster 


0.975 


0.975 


0.975 


0.975 


0.975 


0.975 


C/DOF 


486.05/444 


486.05/444 
















0.72/0.965 


0.71/0.969 


0.84/0.709 


0.86/0.679 



Note: °: Confidence contours are shown in Figure[30] *: Used for the X-ray mass estimate. Galactic absorbtion |W«] = [10 20 cm 2 ] 



global temperature 



[kT] = [keV], solar abundance [Z/Z Q ] = [1], cluster redshift [Zclusi 



■ [1], likelihood function C/DOF, and min. criterion ^J^/ Null hypothesis probability. 




channel anercy (*oV) 



Figure 29. Core spectrum of XMMU J1230.3+1339. The spectrum is fit 
with a single temperature spectral model, extracted from a circle centered 
on the emission peak with a radius of 40"(~315 kpc) for the 0.5 - 7 keV 
band, grouped to include 20 counts per bin. The fit to the minimally binned 
data resulted in an ambient cluster temperature of 5.29 ± 1.02 keV. The Fe 
K-line is visible at 3.3 keV. 



of 2.64 x 10 cm (Dickey & Lockman [T99i3) . and the metal abun- 
dance was fixed at 0.3 solar (Edge & Stewart [l99 if) . Data with en- 
ergies below 0.5 keV and above 7.0 keV were excluded from the 
fit. 

The best fit temperature of the spectrum extracted from 
within this region was T x = 5.29 ± 1.02 keV, with a Cash- 
statistic C = 486.05 for 444 degrees of freedom. Using this spec- 
trum, the unabsorbed luminosity of XMMU J1230.3+ 1339 within 



40"(315 h^ 1 kpc) radius is L x = 1.48 x lO^erg s _1 (0.5-2.0 keV). 
In Figure [29] we show a plot of the spectrum with the best fitting 
model overlaid. The spectrum was grouped to include at least 20 
counts per bin. We performed several tests in order to assess the 
robustness of our results against systematics, by changing back- 
ground regions and by using different statistics. A comparison of 
the Cash-statistic and ^-statistic is shown in Table 4; the best fit 
temperature varies by only ~ 10 %. The impact of a possible cali- 
bration problem of ACIS-S at high-energies has been alleviated by 
using the latest calibration files. 

The resulting temperature, combined with the density profile 
inferred from the ft model fit 

-3,8/2 



PgasO) = P0 



1 + 



and the mass equation derived from hydrostatic equilibrium 



M tot (< r) : 



kT(r)r IS In p <51nT 



G/jnip \ 6 In r 6 In r 



(38) 



(39) 



(where /im p is the mean mass per particle), were then used to esti- 
mate the value of i 2 qo, following Ettori (2000) 



^200 : 



3f3kT 



Gnm„(4/3)np c (z)^200 



1, 



(40) 



where r c is the core radius. With /j = 0.62, we obtain r2oo = 1046 ± 
113 kpc. The total mass is then calculated out to r 2 oo, using the 
equation 

•\2 

(41) 



(r/r c f 



G iim p 1 + (r/r c ) 2 
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Table 5. Mass estimates for XMMU J1230.3+1339. 



Reference 


^200 


M 20 o 




Method 




[kpc] 


[10 14 M Q ] 


[10 14 M Q ] 


(Remark) 


This work 


944 + 397 


2.85 + 3.6 


3.01 + 3.8 


Kinematics 


This work 


1046 ± 113 


3.87 ± 1.3 


3.72 ± 1.2 


X-ray 


This work 


1376 ±218 


8.80 ±4.2 


6.74 ± 3.3 


WL (NFW)° 


This work 


1125 ± 188 


4.56 ±2.3 




combination 


Paper I 


1017 ±68 


4.19 ±0.7 




combination 



°: Note that we use for the ftoo the 50% significance level. 
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Figure 31. Mass Profile of XMMU J1230.3+1339. In black: best fit NFW 
profile from weak lensing; in blue: best fit kinematic SIS profile; in red: 
best fit X-ray /3-model profile. We have added the M900 and f200 values with 
errorbars for each method. 



7 CLUSTER MASS ESTIMATES 

In this section we compare the masses determined with different 
methods, namely from the kinematics of the cluster galaxies, X-ray 
emission of the cluster gas, and weak lensing. 

As seen from Table 5 the mass estimates from the NFW fit, 
kinematics, and X-ray analyses are in good agreement, within the 
error bars. In Figure [3j] we plot the mass estimates against the 
corresponding radius. The dynamical mass estimate is on the low 
side but still consistent with the other measurements due to the low 
statistics. 



8 DISCUSSION 

Here we discuss the results of our optical, weak lensing, kinematic 
and X-ray analyses and compare XMMU J1230.3+ 1339 to other 
high redshift clusters (see e.g. Appendix Table CI). A more gen- 
eral discussion of the global cluster properties (e.g., mass proxies, 
scaling relations, large scale structure) can be found in paper I. 



8.1 Colour-Magnitude Relation 

We observe evidence for a deficit of faint red galaxies, down to a 
magnitude of z AB = 24.5 (50 % completeness), based on our photo- 
z/SED classification (see Figures [751 and [T8l. Similar results have 
been reported for sever al clusters at z = 0.7 - 0.8, for R XJ015 2 
by Tanaka et al. d2005h . for RXJ1716 by Koyama et al. d2007t) . 
and for 18 EDisCS clusters by De Lucia et al. fcoOTt) . and for RD- 
CSJ09 10+54 at z = 1.1 by Tanaka et al. d2008h . However there are 
other cluste rs still under debate, like e.g. MS 1054, where Andreon 
et al. ( 120061) show ed the faint end of the CMR and in contrast to that 
Goto et al. (2005), who measured a deficit at lcr level for the same 
cluster. Some other known clusters seem to have their faint end of 
the CMR in place, like e.g. the 8 c lusters between redshift 0.8 and 
1.27 studied in detail by Mei et al. d2009h and also the "twin" clus- 
ter XMMU J1229+ 0151 at z = 0.975 (Santos et al. l2009h . 

Koyama et al. d2007l) suggested to collect more observations to 
study the deficit of faint red galaxies, since it can be seen only in the 
richest clusters and therefore a larger cluster sample is necessary. 

XMMU J1230.3+1339 is a very rich, massive, high-redshift 
cluster, and this may have some impact on the build-up of the CMR. 
We will explore the build-up of the CMR in clusters at z< 1 and do a 
comparison with semi-analytical modeling of hierachical structure 
and galaxy formation in an upcoming paper. 

8.2 Blue galaxy fraction - Butcher & Oemler effect 

About half of the cluster members has a blue galaxy type, as 
inferred from the photo-z/redshift classification. This is in good 
agr eemen t with the predictio n of D e Lucia et al. d2007t) , Menci et 
al. d2008h and Haines et al. (2009): an increase in the fraction of 
blue galaxies with increasing redshift. 

Note that the photometric redshift estimates for blue galaxies 
are usually worse than those for red galaxies, due to the different 
shape of the SED's (see right panel of Figure 1 1 for details) which 
results in larger uncertainties. This can not be seen in this data set, 
because we only have "red" spectroscopic clusters members (see 
Ilbert et al. 2006 and Brimioulle et al. 2008 for the dependence of 
photometric redshift accuracy on different SED types and magni- 
tude intervals). The larger photometric redshift errors of blue galax- 
ies lead to an under estimation of the blue galaxy fraction due to an 
increased scattering of object into neighboring redshift bins, see 
e.g. the small signal bins z = [0.8 - 0.9] and z = [1.0 - 1.5] in Figure 
M 

8.3 Colour-Colour Diagram 

Focussing on the colour coded cluster members, based on the 
photo-z/SED classification one can easily see that a separation be- 
tween foreground and background for the blue galaxies is hardly 
possible without knowing the redshift. The SSP models, computed 
using PEGASE2 (Fioc & Rocca-Volmerange 1997), nicely repro- 
duce the observed colours of the early type galaxies of the cluster 
and the field. Note that when observing a cluster at redshift ~ 1, the 
flux of the cluster members in the blue bands is very low, especially 
in the U-band, therefore the measurement errors are large. 

E.g., additional NIR observations (J and K band down to J* +4 
and K* + 5) can reveal the existence of an additional intermediate 
age stellar population, and constrain the star formation rate (SFR) 
of all cluster members together with the available data. With these 
data we will also be able to the study slope, scatter and zeropoint 
of the RS in the J - K CMR. To further understand the evolution of 
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cluster elliptical galaxies around z = 1 not only the strength of the 
4000 Angstrom break is essential, but also the measurement of the 
fraction of young and intermediate stellar populations. 

8.4 Mass estimates 

8.4.1 Weak lensing analysis 

Very recently Romano et al. performed the first weak lens- 

ing study using the LBT. They analysed the galaxy cluster Abell 
611 at z ~ 0.3 and demonstrated that the LBT is indeed a suitable 
instrument for conducting weak lensing studies. We observe a weak 
lensing signal for a cluster at a redshift of z ~ 1, the most distant 
cluster so far used in a ground-based weak lensing study. We were 
able to derive a weak lensing mass estimate for the following rea- 
sons. 

1. Deep high-quality multi-band data (U, B, V, r, i, z) have been 
taken under subarcsecond seeing conditions in order to resolve 
faint, small and distant objects (cluster members and background 
galaxies). 

2. We obtained accurate photometric redshifts and SED classifi- 
cations to separate cluster members, foreground and background 
galaxies. This avoids a dilution of the shear signal from cluster 
member and foreground galaxies, which would lead to an addi- 
tional mass uncertainty of ~ 4 % (see Hoekstra l2007l) . 

A conventional background separation using the CMR or colour- 
colour diagram is not possible (see e.g. Medezinski et al. 
since the background galaxies are not distinguishable from clus- 
ter members and foreground galaxies in the colour-magnitude or 
colour-colour space. 

3. The PSF pattern was stable and isotropic over the entire ob- 
servation time, and could be modeled with a 5''' order polynomial. 
The systematic error is almost 1/3 of the signal for angular scales 
considered in the weak lensing analysis. 

We had to reject on average one third of the available data (in each 
band) to produce the stacked images for our scientific exploitation, 
see Section |2~T1 for details. Note that our method (KSB) to mea- 
sure galaxy shapes and e xtract t he shear signal is similar to the one 
applied in Romano et al. J2010f) . 

We used spectroscopic data and the two-point angular-redshift 
correlation analysis to investigate the photometric redshift accu- 
racy. The PSF systematic residual has been tested with the cross- 
correlation function between galaxies and stars. All tests perform 
very well and we conclude that the systematics are minimal. We 
derive a velocity dispersion of cr^Sr = 1308+ 284 km/s by fitting 
a SIS model to the tangential shear profile. From the NFW likeli- 
hood analysis we find a best-fit concentration parameter c = 4.0+j 4 
and a scale radius r s = 344*«2 kpc. The concentration parameter 
is only weakly constrained. Nevertheless, the best-fit value in good 
agreement with the relation of Bullock et al. (2001) and Dolag et 
al. (2004). This results in a total mass of M 2 N ™= (8.80 ± 4.17) x 
1O 14 M . 

Hoekstra J2003I) analysed the effect of the large scale struc- 
ture along the line of sight, which also contributes to the lensing 
signal, and consequently affects the mass estimates. This intro- 
duces a factor ~ 3 for the uncertainties in M 2 oo and a factor ~ 2 
for the concentration parameter c at t he redshift of z = 1 (extrapo- 
lated from Figure 7 in Hoekstra 2003). For de ep HST observations 
of MS 1054-0321 (z=0.89) Jee et al. fc005bl) estimated an addi- 
tional ~14% error for the total cluster mass due to cosmic shear. 
Our analysis is limited by the large statistical error of the measure- 



ment (shape noise), which is the driving uncertainty of our mass 
estimate. 

XMMU 11230.3+1339 is located behind the local Virgo clus- 
ter (z = 0.0038), which covers a sky area of 8 x 8 square degrees, 
while XMMU 11230.3+1339 covers 6x6 square arcminutes. 

The Virgo cluster acts like a mass sheet, with a constant sur- 
face mass density (k = const.) leading to a constant deflection angle 
(a = const.) and no shear (y = 0). Note that the surface mass den- 
sity k is indirect proportional to 2 cn -, . The critical surface ma ss den - 
sity E cn - ( depends on the ratio D J(D d D ds ) (Schneider et al. 1 19961) . 
which leads in the case of the Virgo cluster lens geometry to a high 
value of Eerie Therefore k becomes negligible small for the Virgo 
cluster. 

We observe an elongation in the projected mass distribution 
from the M„ p statistic similar to the distribution of galaxies asso- 
ciated with the main cluster component and the group bullet (see 
Paper I) and the X-ray emission. Currently the statistics and the 
depth are not sufficient to resolve the possible substructure peaks 
associated with the central group merger event. The WL separation 
of the main cluster and the group bullet will require high-resolution 
space-based imaging data. 



8.4.2 Kinematic analysis 

The current analysis is based on only 13 cluster members, resulting 
in a total mass of Mf™= (2.85 + 3.60) x 10 14 Mq. We expect to 
improve this result substantially by including more spectroscopic 
data. At the current level the kinematic mass estimate represents 
only a upper limit of the mass. The analysis can be improv ed fol - 
lowing, e.g., the method described in detail in Biviano et al. d2006h . 



8.4.3 X-ray analysis 

The cluster was recently serendipitously observed with ACIS on- 
board Chandra with a total integration time of ~ 38 ks. In spite of 
the surface brightness dimming due to the high cluster redshift, we 
can see extended X-ray emission out to ~ 0.5 Mpc. This emission 
shows a regular morphology and is clearly elongated in the same 
direction as the distribution of red galaxies and the projected mat- 
ter distribution. Our spectral analysis of the Chandra data yields a 
global temperature of T x = 5.29 ± 1.02 keV, and a total mass of 
Mf m = (3.87 ± 1.26) x 10 14 Mq, in good agreement with Paper 
I. We applied several tests (different background regions, different 
X-ray statistics) to check the robustness of our findings, which we 
conclude are robust against several systematics involved in the X- 
ray spectral analysis. 

Any X-ray study relies on the accuracy achieved for the tem- 
perature and electron density distribution of the ICM and thus mea- 
sure the mass distrib utions with statistical uncertainties below 15% 
(e.g., Vikhlinin et al. 2006). However, the accuracy of X-ray cluster 
mass estimates based on hydrostatic equilibrium is limited by addi- 
tional physical processes in the ICM (e.g., merger) and projection 
effect (e.g., substructure). Recently Zhang e t al. ( 120090 showed that 
for substructure (see Bohringer et al. l2010aO this reflects into ~ 10% 
fluctuations in the temperature and electron number density maps. 
Given the evidence for a cluster bullet and up to six galaxy over- 
densities in the cluster (see Paper I) these possible will introduce 
an additional ~ 20% systematic error in the mass estimate. Never- 
theless, our X-ray analysis is still dominated by the poor statistic of 
the data in hand. 

Note that the X-ray global temperature does not follow the 
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Figure 32. Temperature, 7x, and galaxy velocity dispersion, <x, relation for distant (z > 0.8) galaxy clusters (left panel). Velocity dispersions and temperatures 
are taken from Table Bl. The solid line shows the relation kgTx = fim p cr 2 (e.g., see Bohringer et al. 2010b; Rosati et al. 2002). Red open boxes (black dotted 
triangles) mark galaxy velocity dispersions taken from weak lensing (kinematic) analyses. XMMU J1230. 3+1339 with T\ = 5.3 ± 1 keV and <x = 1308 ± 284 
km/s is represented as large filled red box. Ratio, 0"wL/o"Kin> of the velocity dispersion form the weak lensing and the kinematic analysis against temperature 
Tx (right panel) for the nine distant (z > 0.8) galaxy clusters taken from Table Bl. The mean average ratio is 1.14 ± 0.46. 



T x - cr relation shown in Figure [32] which would lead to a global 
temperature of Tx = 2.59 keV. However together with the other 
high-z cluster it is still consistent with the relation within its larger 
scatter. A si milar re sult has been reported for XMMU J1229+015 1 
(Santos etal. l2009h . 

We observe evidence for a ss 14 % bias in the ratio, 0" WL /o"Ki n , 
of the velocity dispersion form the weak lensing and the kinematic 
analysis against temperature T x for the nine distant (z > 0.8) galaxy 
clusters (right panel of Figure[32]l. For XMMU J1230.3+1339 we 
can infer a ratio trwiVcKin = 1.93 ± 0.42, which is on the upper end 
of all clusters. We want to stress here again the lack of statistics 
for the kinematic analysis, which leads to a large uncertainty in the 
ratio. Within the large error bars, the ratios of the kinematical and 
lensing velocity dispersion seem to be the same. The mean average 
ratio is larger than one (1.14 ±0.46), but it needs a larger data set to 
show whether this difference is significant. If the galaxies with mea- 
sured velocities are still falling into the cluster one would expect a 
velocity dispersion biased high. The velocity dispersion would be 
biased low if non cluster member galaxies from the surrounding 
filaments were considered. Our procedure to eliminate galaxies in 
the determination of the velocity dispersion is however made to 
avoid the latter effect. The ratio of weak lensing to dynamical ve- 
locity dispersion could also be biased high when matter associated 
to groups correlated to the cluster and along the line of sight is ac- 
counted to the cluster in the weak lensing analysis. One caveat in 
estimating the lensing to dynamical velocity dispersion lies in the 
reliability of the errors. For instance some of the clusters have a 
weak lensing velocity dispersion of around 900 km/s with a 1 cr 
error of only 50 km/s, or about 5%. Unless using space data with a 
high density of objects such an accuracy is very difficult to achieve. 
Even more so, if one accounts for the fact that clusters at these 
high redshifts consist of many subclumps, and that the redshift dis- 
tribution of background objects shows cosmic variance, the weak 
lensing measurement will be biased. 



Given the poor X-ray photon statistics we were able to derive 
weak constraints on the metal abundance Z = (0.29 ± 0.24) Z Q of 
the intracluster medium from the detection of the Fe-K line in the 
Chandra spectrum. 

8.4.4 Comparison 

We estimate independently and self consistently values of the clus- 
ter mass within R200. the radius of the sphere within which the clus- 
ter overdensity is 200 times the critical density of the Universe at 
the cluster redshift, from kinematic, X-ray and weak lensing analy- 
ses. The X-ray mass profile, based on the assumptions of isothermal 
gas in hydrostatic equilibrium, is found to be in good agreement 
with the mass density profile derived from the kinematic analysis. 
Given the fact that the weak lensing analysis is performed at the 
feasibility limit of ground based observations and therefore leads 
to larger errors, the weak lensing mass density profile is consistent 
with the other data within the errors. Note that all our individual 
mass estimates (see Table 5) are in good agreement with each other. 
The combined weighted mass estimate M2oo= (4.56 ± 2.3) x 10 14 
Mq for XMMU J1230.3+1339 is in very good agreement with the 
combined mass estimate of Paper I, where a total mass M 2 oo= (4.19 
± 0.7) x 10 14 Mq is estimated. 



9 SUMMARY & CONCLUSIONS 

This paper describes the results of the first LBT multi-colour and 
weak lensing analyses of an X-ray selected galaxy cluster from the 
XMM-Newton Distant Cluster Project (XDCP), the high redshift 
cluster XMMU J1230.3+1339 at z ~ 1. The high-quality LBT 
imaging data combined with the FORS2 spectroscopic data allows 
us to derive accurate galaxy photometry and photometric redshifts. 
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Our main conclusions are: 

• The photometric redshift analysis results in an accuracy of 
Az/(1 + z) = 0.07 (0.04) and an outlier rate rj = 13 (0) %, when 
using all (secure, > 95 % confidence) spectra, respectively. 

• We find an i - z, against z colour-magnitude relation with a 
slope of -0.027 and a zeropoint of 0.93. Our i - z, against z colour- 
magnitude relation shows also evidence for a truncated red se- 
quence at z = 24, while rich low redshift systems exhibit a clear 
sequence down to fainter magnitudes. Similar evidence has been 
observed for other clusters (e.g., EDisCS clusters, see De Lucia et 
al. 2007; RXJ1716.4+6708 at z = 0.81, see Koyama et al. 2007). 

• The cluster clearly follows the Butcher & Oemler effect, with 
a blue galaxy fraction of i b « 50 %. 

• We see an indication of recent star formation in the distribu- 
tion of the early-type galaxies in the U—B against r—z colour-colour 
diagram. 

• The brightest cluster galaxy could be identified spectroscopi- 
cally and from photometry as an early-type galaxy. A displacement 
of this galaxy from the X-ray cluster centre of ~ 100 kpc is ob- 
served. 

• Compared to XMMU J1229+0151 at z = 0.97 (Santos et al. 
2009) we observe a similar i-z against z colour- magnitude relation. 

• Compared to RDCS J0910+5422 at z = 1.10 (Tanaka et 
al. 2008) we observe also several additional clumps, with similar 
colours in the colour-magnitude relation. 

• We detect a weak lensing signal at 3.5 <x significance using 
the aperture mass technique. Furthermore we find a significant tan- 
gential alignment of galaxies out to a distance of ~ 6' or ~ 2 r 2 oo 
from the cluster centre. The S/N peak of the weak lensing analysis 
coincides with the Chandra X-ray centroid. 

• We have demonstrated the feasibility of ground based weak 
lensing studies of massive galaxy clusters out to z ~ 1 . 

• From our X-ray analysis, we obtain a global temperature of T x 
= 5.29 keV and a luminosity of L x [0.5-2.0keV] = 1.48 x 10 44 erg 
s~\ within an aperture of 40", in good agreement with the results 
presented in Paper I. 

• Mass estimates inferred from the kinematic, weak lensing and 
X-ray analyses are consistent within their uncertainties. We derive a 
total mass of M 200 =(3.9 ± 1.3) x 10 14 M©, r 200 =1046 ± 113 kpc for 
the X-ray data, which has the smallest uncertainty of the three mass 
estimates. From the weak lensing NFW likelihood analysis we find 
a concentration parameter c = 4.0+ 2 4 and a scale radius of r s = 
344^1) kpc, which is consistent with expectation from simulations 
(Bullock et al. 2001; Dolag et al. 2004). 

• We see evidence for a ss 14 % bias in the ratio, (Twh/cKm, of 
the velocity dispersion form the weak lensing and the kinematic 
analysis against temperature T x for distant (z > 0.8) galaxy clus- 
ters. 

We have been awarded 75 ks ACIS-S onboard Chandra and 
an extensive spectroscopic follow-up on ~ 100 photo-z selected 
cluster-members, using F0RS2 to complement our data set. Both 
will be observed in 2010/1 1. With this additional data we will im- 
prove the characterization of this remarkable rich, X-ray luminous 
galaxy cluster. 
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Figure Al. Absorbed ICM model spectrum of the bremsstrahlung con- 
tinuum and line emission for a plasma at a redshift of z = 0.975 with 
Z/Z G = 0.3 and temperature of T = 5.29 keV observed through a galac- 
tic hydrogen column density of Nh = 2.64 X 10 22 cirr 2 . The redshifted Fe 
K-line (6.7 keV at rest-frame) is visible at 3.3 keV. 
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APPENDIX A: ICM MODEL SPECTRA 

An absorbed ICM model spectrum of the bremsstrahlung con- 
tinuum and line emission for a plasma, representing the galaxy 
cluster is shwon in Figure Al. 



APPENDIX B: PSF ANISOTROPY PATTERN 

In Figure IBTI we plot the measured PSF patterns of all 6 bands, 
represented by the unsaturated stars imaged on the field at high SN. 
The measured PSF widths of stars on the stacks are summarised 
in Table 1. The anisotropy pattern of LBC is similar in terms of 
orientation and size to Suprime-Cam shown in Figure 3 of Okabe 
& Umetsu (2008). 



APPENDIX C: HIGH-Z CLUSTER 

In Table C 1 we give an overview of recently studied high-z galaxy 
clusters. So far only 10 of them were studied with weak lensing. 
Note that only XMMU J1230.3+1339 (z~l), C1G J1604+4304 
(z-0.9), RX J1716.6+6708 and C1G 1137.5+6625 (z~0.8) have a 
weak lensing detection based on ground based data. 
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Figure Bl. PSF anisotropy pattern of the image stacks in the 6 (U, B, V, r, i, z) filters. The sticks mark the observed ellipticities of high signal to noise but 
unsaturated stars in the image stacks. The direction of the sticks coincides with the major axis of the PSF, the length is indicating the degree of anisotropy. An 
ellipticity of 10% is shown in the upper left corner of each panel. The definition of ellipticity is given in Eq. 4. The measured PSF widths of stars on the stacks 
are summarised in Table 1. 



Table CI. List of high-z galaxy clusters and their basic optical/X-ray/weak lensing properties: (1) sequence number, (2) right ascension, (3) declination, (4) number of spectroscopic objects used to estimate the 
cluster redshift (6) and the line-of-sight velocity dispersion (7), (8) bolometric X-ray luminosity, (9) ICM temperature, (10) velocity dispersion from a weak lensing analysis, more information on peculiar individual 
clusters is provided in the associated reference (11). 



Id 

(1) 


Cluster name 

(2) 


R.A. 
(J2000) 

(3) 


Dec. 
(J2000) 
(4) 


N s«' 
(5) 


Redhift 

(6) 


of" 
[km/s] 
(7) 


„, L * D ' , 
[lO^ergs -1 ] 

(8) 


T v 
[keV] 

(9) 


[km/s] 
(10) 


Ref. 
(11) 


1 


SXDF-XCLJ0218-0510 


02 , T8'"21.3 1 


-05°10'27" 


9 


1.6230 + 0.001 


537 ± 


213 


0.34 ±0.16 


1.7 ±0.3 


... ± ... 


Tanaka et al. (2010); Papovich et al. (2010) 


2 


XMMXCS 12215.9-1738 


22*15 m 58.5 s 


-17°38'03" 


31 


1.457 ±0.001 


720 ± 


110 


2.92 ± 0.35 


4.1 ±0.9 


... ± ... 


Hilton et al. (2010); Stanford et al. (2006) 


3 


XMMU 12235.3-2557 


22' ! 35 m 20.8 I 


-25°57'40" 


30 


1.390+ ... 


802 ± 


77 


8.5 ± 0.4 


8.6 ± 1.3 


1 145 ± 70 


Jee et al. (2009a); Rosati et al. (2009); 
Mullis et al (90051 

IVlLllllj CL til. \ Z_V'V— ' ) 


4a 
4b 


CIO 10848+4453 1 vnx W 
C1G J0849+4452, Lynx E 


08 ,! 48 m 34 7 s 
08*48 m 56.3 s 


+44°53'42" 
+44° 52' 16" 


9 
11 


1.273 ± 
1.261 ± ... 


650 ± 
720 ± 


170 
140 


1.5 ± 0.8 
2.1 ±0.5 


1.7 + 0.7 

3.8 ± 1.3 


762 ± 133 
740 ± 134 


Mei pr al c 900QV 
Jee et al. (2006) 


5 


RDCS J 1252.9-2927 


12' ! 52 m 48.0 I 


-29° 27' 00" 


38 


1.2370 ± 0.0004 


747 + 


84 


6.6 ± 1.1 


6.0 ±0.7 


1185 ± 120^ 


Demarco et al. (2007); Lombardi et al. (2005); 
Rosati et al. (2004) 


6 


RDCS J09 10+5422 


09 ,! icqo-o 1 


-54°22'00" 


161 


1.1005 ±0.0016 


716 ± 


141 


2.9 ± 0.4 


2.9 ± 0.4 


... ± ... 


Maughan et al. (2008a); Tanaka et al. (2008) 


7 


XI SSI072404 1-041330 


02''24 m 04 I s 

\Jj— AST V/^T. 1 


-04° 13'30" 


5 


1.050 ± 


... ± 




4.5 ±0.1 


4.3 ±1.1 




Mmiohan pt al f?008hV Paranrf pt al (90071 

IVluU £fll£lll CI ell . ^UOGU), L uCuU-U CI dl . ^ iL\J\J 1 J 


X 




02 /! ?l m 42 s 


-03°21'47" 


5 


1.016 ± 


710 ± 


150 


+ 


+ 


+ 


AnHrenn pt al f900R1 

^illLllCUll CL dl. V ) 


9 


XMMU 11229.4+0151 


12*29 m 29.2 s 


+01°51'26" 


27 


0.975 ± 0.001 


683 + 


62 


3.0 ± ... 


6.4 ± 0.7 


... ± ... 


Santos et al. (2009) 


10 


XMMU 11230.3+1339 


12*3Q"'16.0' 


+ 13°39'04" 


13 


0.9745 ± 0.002 


657 ± 


277 


6.5 ± 0.68 


5.29 ± 1.02 


1308 ± 284 


This work & Paper 1 


11a 
lib 


C1G J 1604+4304 
C1GJ1604+4314 


16*04 m 23.7 s 
16 /, 04 m 25.8 I 


+43°04'52" 
+43°14'23" 


67 
62 


0.9001 ± ... 
0.8652 ± ... 


962 ± 
719 ± 


141 
71 


1.576 ± ... 
1.164 ± ... 


3.5 ± 1.82 

1.6 ±0.65 


1004 ± 199 
... ± ... 


Kocevski et al. (2009); Margoniner al. (2005); 
Gal et al. (2004) 


12 


C1G J 1226+3332 


12' ! 26 m 58.r 


+33°32'47" 


6 


0.8877 ± 0.004 


1650 ± 


930 


5.12 ± 0.12 


10.4 ± 0.6 


880 ± 50 


Jee et al. (2009b); Maughan et al. (2007); 
Ebeling et al. (2001) 


13 


RX 11257.2+4738 


12*57«12.2 S 


+47°38'07" 


18 


0.866+ ... 


600 i 




... ± ... 


3.6 ± 2.9 


... ± ... 


Ulmer et al. (2009) 


14 


C1G 1054-0321 


10*57 m Q0.2 s 


-03°37'27" 


129 


0.826+ ... 


1156: 


t 82 


3.1 ±0.1 


8.9 ± 1.0 


1150 ± 51 


Tran et al.(2007); Jee et al. (2005b); 
Hoekstra et al. (2000) 


15 


C1G JO 152- 1357 


01*52 m 41.0 s 


-13°57'45" 


102 


0.837 ± 0.001 


1632 ± 


115 


... ± ... 


6.46 ± 1.74 


903 ± 57 


Jee et al. (2005a); Demarco et al. (2005); 
Maughan et al. (2006) 


16 


RXJ1716.6+6708 


17 ,1 16 m 49.6 I 


+67°08'30" 


37 


0.809 ±0.0051 


1522 ± 


215 


1.50 ±0.4 


5.7 ± 1.4 


1030 ± ... 


Clowe et al. (2000); Gioia et al. (1999) 


17 


C1G 1137.5+6625 


ll' ! 40'"23.3 I 


-66°09'01" 


22 


0.784 ± 0.0003 


884 ± 


185 


1.90 ±0.4 


5.7 ±2.1 


1190 ± ... 


Clowe et al. (1998); Luppino et al. (1995); 
Donahue et al. (1999) 



For a non-singular isothermal sphere (NSIS) with free coordinates for its center. 



